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The dragonfly genus Celithemis consists of 8 spp., some of them common and 
brightly colored, that are confined largely to eastern North America. Several spp. have 
been used in behavioral, ecological, and morphological studies, but their intrageneric 
phylogeny is unclear. In this paper is provided a phylogeny based on morphology and 
on data from mitochondrial and nuclear DNA sequences of multiple individuals of 
each species. The genus appears to be monophyletic, with one nested species pair (C. 
amanda + C. martha) receiving strong bootstrap support by both parsimony or maxi¬ 
mum-likelihood criteria as well as high Bayesian posterior probability. A second group 
(C. bertha, C. elisa, C. ornata and C. fasciata ) is well-supported in Bayesian analysis 
but only weakly by parsimony and maximum-likelihood bootstrap values. C. verna 
and C eponina are probably basal to both these groups, but their relationship to each 
other is unclear. All individuals assigned to a species recognized on morphological 
-grounds were recovered as monophyletic. The problematic taxa, C. monomalaena and 
C. bertha leonora , are shown definitively to be synonyms of C. fasciata and C. bertha , 
respectively. 


INTRODUCTION 

Aside from the inherent interest of their evolutionary history, Odonata have 
been important in research on freshwater systems as well as behavioral ecology, 
sexual selection, and insect physiology and functional morphology. The genus 
Celithemis , in particular, has been used in studies of functional morphology of 
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the penis (MILLER, 1981), reproductive behavior (MILLER, 1982), competition 
(COWLEY & JOHNSON, 1992), energy balance (BENKE & BENKE, 1975, 
BENKE, 1976), and predation pressure (GRESENS et al., 1982). Understand¬ 
ing relationships among taxa is essential for comparative assessment of adaptive 
similarities and differences in behavioral ecology. Thus establishment of a reliable 
phylogeny is a prerequisite for any deep understanding of selection and character 
evolution in Celithemis , as in any group of organisms. 

HAGEN (1861) established the genus Celithemis to include Libellula eponina 
Drury, Celithemis superbum sp. n. and the fossil species, C. cellulosa sp. n. Ce¬ 
lithemis superbum was moved to the monotypic genus, Pseudoleon Kirby, 1881, 
and C cellulosa is clearly distant from our current concept of Celithemis (see 
GARRISON et al., 2006, for a full diagnosis) by virtue of having the midrib of 
the anal loop strongly bent, and having many more cells in the forewing triangle 
and supratriangle and many more antenodal crossveins than in extant species. 
Thus these two species are not considered further. 

Celithemis has not been revised since WILLIAMSON (1922), and no mod¬ 
em analysis of the species boundaries or relationships has been attempted. The 
named taxa, C. monomalaena Williamson, 1910, and C. bertha leonora Westfall, 
1952, have recently been considered synonyms of C. fasciata Kirby, 1889, and C. 
b. bertha Williamson, 1922, respectively, but unequivocal data supporting syn¬ 
onymy have not been published. KORNELL (2003) investigated phylogeny and 
species distinctions using molecular and morphological data, but several gaps 
and inconsistencies in her data have encouraged this re-examination. 

Our principal objective here is to derive a well supported phylogeny resolving 
the relationships among species of Celithemis . Since morphological features lead 
to little resolution, we combine this with sequence data from the mitochondrial 
protein coding gene cytochrome oxidase I (COI) and from the D2 region of the 
nuclear 28S rDNA gene. We also provide data supporting the synonymy of C. 
monomalaena and C bertha leonora . 

MATERIAL AND METHODS 

MOLECULAR TECHNIQUES — Twenty-eight specimens for molecular analysis, including rep¬ 
resentatives of all putative species and subspecies of Celithemis and three outgroup species, Cordulia 
shurtleffi, Libellula pulchella, and Leucorrhinia glacial is, were field collected or taken from the collec¬ 
tions of MLM and Rutgers University. Tissue was obtained by excising thoracic muscle or removing 
a leg, which was then placed in an Eppendorf tube. Extraction procedures followed manufacturer’s 
recommendations (QIAGEN, DNeasy Tissue handbook). 

After the DNA was extracted, polymerase chain reaction (PCR) was used for amplification of a 
COI fragment from the mtDNA and the D2 region of the large subunit (28S) nuclear ribosomal 
coding region. We used the C1731F and CI-14R primers of Artiss et a/.(2001) for COI amplification 
and the D20DUP and D2DnB primers (5’TGCTTGAGAGTGCAGCCCAA3’ and 5’CCTT- 
GGTCCGTGTTTCAAGAC3’). We amplified the DNA with 1 pi each of up and down primer, 
1 to 5 pi of template DNA, 12.5 pi Master Mix solution, and 5.5 to 9.5 pi ddH,0. 



Phylogeny of Celithemis 


103 


Amplified DNA was separated on a 1.5% low melting agarose gel (FMC Bioproducts), then pu¬ 
rified using QIAquick PCR Purification Kit Protocol (QIAGEN) according to the manufacturers 
instructions. The purified product was sequenced in both directions using BigDye Terminators Cy¬ 
cle Sequencing Kit (Applied Biosystems) and an ABI 3100 capillary sequencer. COI sequences were 
aligned using Clustal (THOMPSON, 1997) then, to assure homology, translated to the correspond¬ 
ing protein sequence in MacClade (M ADDISON & MADDISON, 1992). D2 sequence alignments 
were made using Clustal, and the resulting files were then aligned manually in Microsoft Word using 
the structural methods described in KJER et al. (1994), KJER (1995), KJER et al., (2007) and sec¬ 
ondary structure models based on GUTELL et al. (1993). 

D2 sequences only were provided for Leucorrhinia orientalis (specimen obtained from Dr Elena 
Malikova) and a Leucorrhinia sp sequence obtained from Gen Bank (Acc. # AY859583)). COI se¬ 
quences were scored as missing data for these two species. 

Thirty-six morphological characters, including 3 larval characters, were scored for all putative taxa 
of Celithemis and for the 3 outgroup species. All adult characters were observed directly, mostly us¬ 
ing a Wild™ stereomicroscope. Characters of the secondary genitalia were examined with an Hitachi 
S510 scanning electron microscope after sputter coating with gold-palladium. Data on larval char¬ 
acters were taken from NEEDHAM et al. (2000). 

PHYLOGENETIC ANALYSIS — Each data partition was analyzed individually and in com¬ 
bination. The partitioned analyses were not seriously incongruent, although they differed in their 
level of resolution. Morphology and D2 gave poor resolution, while the COI tree was well resolved. 
The analyses of the combined data are presented hereafter. A maximum parsimony heuristic search 
of 1,000 random replicates was conducted using PAUP- 4.0b 10 (SWOFFORD, 1999) with TBR 
branch swapping. To assess branch support, 400 bootstrap pseudoreplicates (FELSENSTEIN, 1985) 
were performed using 10 random addition searches per pseudoreplicate. Modeltest 3.06 (POSADA 
& CRANDALL, 1998) was used to obtain an evolutionary model for each molecular data set in the 
Bayesian and maximum likelihood analyses; in each case the GTR+I+ T was estimated to be the 
best model for molecular data, with mk model used for morphology. Using the program MrBayes 
(HUELSENBECK, 2000), Bayesian analysis (4 MCMC chains: 1 cold, 3 hot) was run on the com¬ 
bined data for 2,873,500 generations; the first 100,000 were discarded as bum-in. 

MORPHOMETRIC ANALYSIS - WILLIAMSON (1910) distinguished C. monomalaena from 
C. fasciata based on size, wing coloration, and several venational characters. Numerous authors (e.g. 
NEEDHAM et al., 2000) have doubted the validity of C. monomalaena , and CARLE (1982) indi¬ 
cated that the supposed differences were either clinal or showed discordant variation. Nonetheless, no 
quantitative evaluation of the supposed distinctions has been published. We measured to the nearest 
0.5 mm or counted the relevant characters (see Results) in 59-66 specimens from the National Mu¬ 
seum of Natural History (NMNH), the Florida State Collection of Arthropods (FSCA), and the 
collection of MLM, using a Wild stereomicroscope when necessary. Data were plotted as a function 
of latitude and evaluated by linear regression to assess geographic variation of each character. 

RESULTS 

Figure 1 presents a strict consensus of the four most parsimonious trees (length 
1065). Celithemis is recovered as monophyletic with good support. All named taxa 
are recovered except that the two C. amanda are in a polytomy with C martha. 
Note also that C. fasciata and C. monomelaena are mutually polyphyletic and 
C. b. bertha is paraphyletic relative to C. bertha leonora. C. verna and C eponina 
form a polytomy with an apparent clade that includes all other species (but this 
clade is supported with very low bootstrap support). The strongly supported C 
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amanda and C. martha clade is apparently basal in the pectinate arrangement of 
the remaining species. 

Figure 2 is the majority rule consensus tree obtained by Bayesian analysis. Ce- 
lithemis verna is sister to the remaining species, followed successively by C. eponina 
and the C. amanda + C. martha clade and finally a polytomy comprising C. ber¬ 
tha, C. ornata + C. elisa, and C. fas data. Posterior probabilities provide strong 
support for the monophyly of Celithemis and for that of most of the species ex- 



■ Cordulia shurtleffii 

• Libellula pulchella 

• Leuchorrhinia glacialis 

• Leucorrhinia orientalist 

■ Leucorrhinia sp . 2 

■ Celithemis verna, MD 

■ Celithemis verna, MD 

■ Celithemis verna, NJ 

■ Celithemis verna, NJ 

■ Celithemis eponina, MD 

■ Celithemis eponina, NJ 

■ Celithemis eponina, TX 

• Celithemis eponina, FL 

■ Celithemis eponina, NJ 

1 Celithemis amanda, SC 
Celithemis amanda, SC 

1 Celithemis martha, NJ 

' Celithemis martha, NJ 

’ Celithemis martha, NJ 
Celithemis martha, NJ 

■ Celithemis ornata, FL 
Celithemis ornata, FL 
Celithemis bertha bertha, FL 
Celithemis bertha bertha, SC 
Celithemis bertha "leonora", SC 
Celithemis elisa, SC 
Celithemis elisa, TN 
Celithemis elisa, KY 
Celithemis fasciata, SC 
Celithemis "monomelaena", NJ 
Celithemis fasciata, SC 
Celithemis fasciata, AL 
Celithemis " monomelaena", NJ 
Celithemis " monomelaena ", NJ 
Celithemis fasciata, TX 
Celithemis fasciata, AL 


Fig. 1. Strict consensus of four most parsimonious trees of length 1065 based on combined data 
from COI, D2 and morphology. Numbers above the branches indicate bootstrap values > 50%. The 
names in quotations, “leonora” and “monomalaena”, are regarded as synonyms of C. bertha and C. 
fasciata , respectively, but are retained in the tree to indicate specimens that exhibit characteristics of 
those supposed taxa. Superscripts indicate specimens for which sequence data were obtained from 
E. Malikova (1) or from GenBank (2). 
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cept C. amanda, which is, however, recovered as monophyletic. C. fasciata and 
C. monomelaena are again mutually polyphyletic, and C. b. leonora forms a clade 
with C. b. bertha . Species groupings are at least moderately well supported ex¬ 
cept for the terminal polytomy and within the latter, the possible clade including 
C. ornata and C elisa. It is worth noting that Bayesian analysis of the COI data 
alone (which was the most informative data partition used) recovered C. verna, 
C. eponina , and C. amanda + C martha in a basal polytomy but gave fairly strong 
support (93% posterior probability) to the C. ornata + C. elisa clade. 

Morphometric analyses of specimens of the C. fasciata!C. monomalaena clade 
show gradual latitudinal dines in hindwing length, number of cells between veins 
A 2 and A 3 in the hindwing, number of hindwing marginal cells basal to the anal 



Cordulia shurtleffii 
Libelulla pulchella 
Leucorrhinia sp.2 
Leucorrhinia orientals' 
Leucorrhinia gladalis 
Celithemis verna, MD 
Celithemis verna, MD 
Celithemis verna, NJ 
Celithemis verna, NJ 
Celtihemis eponina, MD 
Celithemis eponina, NJ 
Celithemis eponina, FL 
Celithemis eponina, TX 
Celithemis eponina. NJ 
Celtihemis amanda, SC 
Celtihemis amanda. SC 
Celtihemis martha. NJ 
Celtihemis martha, NJ 
Celtihemis martha. NJ 
Celtihemis martha, NJ 
Celtihemis bertha bertha. SC 
Celtihemis bertha 'leonora", SC 
Celtihemis bertha bertha, FL 
Celithemis ornata, FL 
Celithemis ornata, FL 
Celithemis elisa. SC 
Celithemis elisa. TN 
Celithemis elisa. KY 
Celithemis ‘monomelaena", NJ 
Celithemis ‘monomelaena", NJ 
Celithemis fasciata, TX 
Celithemis fasciata, AL 
Celithemis fasciata, SC 
Celithemis fasciata, SC 
Celithemis fasciata, AL 
Celithemis ' monomelaena ' NJ 


Fig. 2. Majority rule consensus of 11,096 trees obtained for combined COI, D2, and morphology data 
from Bayesian analysis. Posterior probabilities > 50% are shown above the corresponding branches. 
See Fig. 1 for explanation of quotation marks and superscripts. 
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loop (Fig. 3A), and extent of the forewing nodal spot; hindwing length also tends 
to increase clinally from east to west (p = 0.09, not shown). Characters claimed 
by WILLIAMSON (1910) to distinguish the taxa do not vary concordantly from 
north to south (Fig. 3B), although specimens from north of 38 N do tend on aver¬ 
age to have fewer forewing triangle and posttrigonal cells and fewer basal paranal 
cells (see NEEDHAM et al., 2000). The proportion of specimens with yellow at 
the base of the hindwing, supposedly characteristic of nominate C. fasciata , was 
not quantified because of variation in its intensity and extent, but, although this 
color pattern was clearly more frequent in southern specimens, individuals with 
completely clear hindwings occurred as far south as 32.3° N (GA, Crawford Co.) 
and with distinct yellow basal areas as far a 42.2° N (NY, Chataqua Co.). 

DISCUSSION 

All 8 generally acknowledged species appear to be valid, monophyletic taxa, 
with the possible exception of C. martha , which could be a junior synonym of C. 
amanda ; this pair is quite similar in general appearance and morphology except 
for the strong darkening of mature males of C. martha. We do not advocate a 
taxonomic change, however, pending a more thorough study of these two nomi¬ 
nal taxa. 

We suggest that the putative taxa C. monomalaena and C bertha leonora be rel¬ 
egated to synonymy with C. fasciata and C. bertha , respectively. In the first case, 
individuals identified as C. monomalaena on morphological grounds are mingled 
haphazardly in a clade with typical C. fasciata (Figs 1-2). Furthermore, Figure 3 
provides no morphological basis for distinguishing these taxa even as subspecies. 
In the case of C. b. leonora and C. h bertha , recognition of “leonora” leaves the 
nominotypical form paraphyletic (Fig. 1). As already described by WESTFALL 
(1952), the subterminal dark wing spot distinguishing “leonora” varies from 
about 0.5 cm in diameter to vanishingly small, so that the supposed subspecies 
is doubtfully diagnosable. Moreover, the “leonora” form is not geographically 
distinct by virtue of being restricted to the Florida panhandle but also occurs at 
least in South Carolina, from which our sequenced specimen originated. Thus it 
is geographically broadly coextensive with C. b. bertha and hence does not fit the 
usual definition of a subspecies (e.g. MAYR, 1963). 

No formal morphology-based phylogeny has been proposed nor, to our knowl¬ 
edge, has any other discussion of relationships within the genus appeared since 
WILLIAMSON’S (1922) very brief treatment. One might assume a close rela¬ 
tionship among C. elisa, C. eponina , and C. fasciata based on wing coloration, as 
Williamson did, or among C. eponina, C. fasciata, C. ornata , and C. verna based 
on larval features (NEEDHAM et al., 2000), especially the unusual eye morphol¬ 
ogy of C. verna, C. fasciata, and C. ornata , but neither of these groupings is sup¬ 
ported here. Rather, elaborately patterned wings and conical larval eyes and long 
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abdominal spines are likely to have evolved independently multiple times. A strict 
consensus tree (not shown) based on morphology alone and analyzed using the 
maximum parsimony criterion showed little resolution within Celithemis except 
that it grouped C. amanda and C. martha, based on thoracic color and the shape 
of the cornu of the penis, and C. elisa, C. eponina , and C fasciata based on wing 
coloration. 

Our hypothesis is not without weakness. The parsimony tree (Fig. 1), although 
well-resolved, provides minimal support for interspecific relationships. The Baye- 




Fig. 3: (A) Relation of cell counts and hindwing length, suggested by WILLIAMSON (1910) to dif¬ 
ferentiate G fasciata and C. monomelaena , as a function of latitude. Lines represent linear regres¬ 
sion of the respective measurements on latitude. Each quantity declined significantly (p < 0.001) but 
smoothly with latitude, indicating a N-S cline; measurements of the transverse extent of the forew¬ 
ing nodal spot gave similar results (not shown; p < 0.01). — (B) Proportion of cell counts, suggest¬ 
ed by Williamson to characterize the supposed northern species, C. monomelaena , as a function of 
latitudinal range. In each case the highest proportion occurs at the highest range but with secondary 
peaks at the lowest range and with minima discordant, at 30-34°N (marginal cells and postrigonal 
cells) or 34-38 °N (cells in triangle). 
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sian tree (Fig. 2), however, gives high posterior probabilities for most branches, 
and the two trees are concordant in most respects. The relative positions of C. 
eponina and the C. amanda + C. martha clade, and the relationships among C 
bertha , C. ornata, C elisa and C. fasciata are uncertain, but in other respects we 
think the phylogenies have numerous common features that can be accepted with 
some confidence. 

We hope this study will encourage additional comparative investigations of this 
interesting genus. To mention only two intriguing possibilities: (1) since extensive 
wing maculation at or beyond the nodus appears to have evolved independently 
two or three times, studies of the effect of this character on intra- or interspecific 
signaling could be productive; (2) similarly, extremely conical larval eyes seem to 
have evolved independently two or three times, so the adaptive advantages and 
physiological consequences of this unusual morphology are open to investigation. 
Although the genus is a small one, it may well provide numerous other opportu¬ 
nities for comparisons of behavioral, morphological, and ecological properties. 
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